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The imbalanced antiferromagnet in an optical lattice 
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We study the rich properties of the imbalanced antiferromagnet in an optical lattice. We present 
its phase diagram, discuss spin waves and explore the emergence of topological excitations in two 
dimensions, known as merons, which are responsible for a Kosterlitz-Thouless transition that has 
never unambiguously been observed. 
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Introduction. — Owing to their exquisite experimen- 
tal tunability, ultracold atomic gases have become a test 
bed for many paradigmatic ideas in quantum many-body 
physics. In particular, ultracold atoms confined in an op- 
tical lattice, i.e., a periodic potential for neutral atoms 
created by orthogonal retroreflected laser beams, are ac- 
curately described by the Hubbard model. The nar- 
row energy bands resulting from the periodic potential 
quench the kinetic energy of the atoms with respect 
to their interaction energy, enabling the exploration of 
strongly correlated phases that play a significant role 
in condensed-matter physics. The single-band Hubbard 
model is realised by cold atoms when the lattice potential 
is sufficiently strong that only the lowest-energy band is 
populated 1] . At a filling corresponding to one particle 
per lattice site, this model has a Mott-insulator phase 
for repulsive interactions, and in an atomic Bose gas 
the theoretically predicted superfluid-to-Mott-insulator 
phase transition ^ has indeed been observed experi- 
mentally For bosons, one commonly refers to this 
model as the Bose-Hubbard model. The Fermi-Hubbard 
model, referred to simply as the Hubbard model, has also 
been realised experimentally 0] , and the fermionic Mott- 
insulator has only recently been seen 0] . 

Aside from being rich in novel physics in itself, this 
latter system may also shed new light on the poorly 
understood phenomenon of high-temperature supercon- 
ductivity in the cuprates. There, the Hubbard model is 
thought to describe electrons in the periodic ion-lattice 
potential of the copper-oxygen planes which are believed 
to undergo a quantum phase transition to a d-wave su- 
perconducting state as the filling fraction is reduced by 
doping Q. In this case, electrons on the same lattice site 
repel each other due to the Coulomb interaction. We 
therefore consider here the repulsive Hubbard model, al- 
though the attractive Hubbard model can also be realised 
and possesses interesting BEC-BCS crossover physics Q. 

While achieving the antiferromagnetic ground state of 
the Hubbard model in an atomic Fermi gas is currently a 
major experimental goal, another prominent topic that is 
the focus of a number of recent ground-breaking experi- 
ments is imbalance in two-component Fermi gases 0, [l^ ■ 
There, the population of each spin species can be con- 
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FIG. 1: Phase diagram of the imbalanced antiferromagnet in 
three dimensions. The solid curve is the critical temperature 
ksT / J for the canted antiferromagnet (C). Above the criti- 
cal temperature the system is paramagnetic (PM). Here, J is 
the superexchange coupling of the antiferromagnetic Heisen- 
berg model. Also shown is the critical temperature for the 
corresponding Ising model (dashed line). In the absence of 
imbalance the system is in a pure antiferromagnetic or Neel 
(N) state. 



trolled permitting the exploration of uncharted imbal- 
anced phases of great interest in condensed-matter, nu- 
clear, high-energy, and astroparticle physics. In this Let- 
ter, we explore the bridge between these two exciting 
directions in the physics of ultracold atomic gases. 

Phase diagram — In the Mott-insulator phase only 
spin degrees of freedom remain and its low-lying exci- 
tations are then effectively described by the Heisenberg 
model. The three-dimensional phase diagram calculated 
within mean-field theory for this model is shown in Fig.[T] 
for a cubic lattice which is experimentally most relevant. 
When there are equal proportions of each spin species, 
the Mott insulator becomes a pure antiferromagnetic or 
Neel state below a certain critical temperature and is 
solely characterised by a nonzero expectation value of the 
Neel order parameter vector n, also known as the stag- 
gered magnetisation. Changing the proportion of each 
spin species gives rise to a nonzero average magnetisa- 
tion m = (0,0, TO^), where = S[N^ -Ni)/{N^+Ni) 
is equal to the average spin per site of the system. Here, 
S = 1/2 is the atomic pseudospin for the two-component 
Fermi mixture of interest, and iV|,x denotes the number 
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FIG. 2: The dispersion relation for the spin waves. Here, 
k = fc(l,l,l) is the wave vector and d is the lattice spac- 
ing. For zero magnetisation (solid line), the dispersion is the 
usual doubly degenerate gapless antiferromagnetic dispersion. 
Imbalance splits the degeneracy and one of the dispersions 
(dashed line) becomes ferromagnetic and acquires a gap equal 
to 2Jzmz where z is the number of nearest neighbours, while 
the other remains antiferromagnetic (dotted line). The im- 
balanced dispersions shown here are for — 0.2. 



of particles of each spin species. In this imbalanced case, 
the system has a canted antiferromagnetic phase at low 
temperatures where the Neel vector is always found to 
be perpendicular to the magnetisation 11]. Imbalance 
therefore breaks the rotational symmetry of the system 
and defines an easy plane for the Neel vector that is per- 
pendicular to the average magnetisation. 

In both the Neel and canted phases the low-energy 
excitations are spin-density waves called magnons. By 
linearising the Heisenberg equations of motion of the on- 
site spin operators, we find that for the balanced situation 
the magnons have the usual doubly degenerate antifer- 
romagnetic dispersion that is linear for small momenta, 
as shown in Fig. [21 However, in the imbalanced case the 
nonzero average magnetisation lifts this degeneracy. As 
a result, one of the dispersions becomes gapped corre- 
sponding to Larmor precession of the spins around the 
effective magnetic field generated by the nonzero aver- 
age magnetisation. The gapped dispersion is quadratic 
at small momenta signalling the appearance of magnons 
with a ferromagnetic character. 

In experiments the atomic gas is always confined in 
a trap. However, including a smooth harmonic trap- 
ping potential does not lead to inhomogcncitics in the 
Mott-insulator state but places a limit on the total num- 
ber of particles beyond which the Mott insulator is de- 
stroyed [l^ . As a result, the excitation spectrum shown 
in Fig. [5] becomes discretised in the trap due to finite-size 
effects. In the Mott-insulator phase, the system in a trap 
is thus essentially hornogeneous and we may reasonably 
neglect the trap here |13j . 

Merons — The long-wavelength dynamics obtained 
above can be summarised by the following non-linear 
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Here z = 2D is the number of nearest neighbours in the 
D-dimensional hypercubic lattice, d is the lattice spacing, 
m = (0, 0, TTiz) is the average magnetisation per site, and 
n(x, t) the local staggered magnetisation at position x 
and time t with fixed length n. The equilibrium value 
of the staggered magnetisation is found from minimising 
the Landau free energy, given by 

i^[n(x),m]= I ^|:^[Vn(x)]V/[n(x),m]|, 

(2) 

where /[n(x),m] is the on-site free energy which, as we 
have seen, breaks the rotational symmetry of the system 
in the presence of imbalance. 

The equation of motion that follows from the non- 
linear sigma model exhibits long-wavelength spin waves 
discussed above but also admits a class of topologically 
stable excitations. In the three dimensional case, these 
are coreless vortices with a spin component near the lo- 
cation of the vortex line that is out of the easy plane. 
In the two-dimensional case the topological excitations, 
known as merons are particularly interesting. 

Far away from the core of the meron, n(x) lies in the 
easy plane and forms a vortex with vorticity of ±1, while 
in the core region n(x) smoothly rotates either up or 
down out of the easy plane. This spin texture is char- 
acterised by two topological invariants, namely, the vor- 
ticity and the Pontryagin index which is equal to ±1/2 
depending on the polarity of the core spin [17| . The Pon- 
tryagin index does not play an important role for our 
purposes since, due to the symmetry of the action un- 
der the reflection n —^ — n, merons differing only in their 
Pontryagin index have the same energy. Moreover, the 
interaction between two merons with a large separation is 
insensitive to the precise spin configuration in the core. 
The energy of a single meron diverges logarithmically 
with the system area A as ( Jn^7r/2) ln(A/7rA^). Here, 
A = dn^JirJ l2Fm is the characteristic size of the core, 
plotted in Fig. [S^b), with Fm the integrated on-site free 
energy required for the formation of a meron spin texture 
of size \ = d. Thus, for large system sizes of interest here, 
it is impossible to thermally excite a single meron below 
the critical temperature. This behaviour is analogous to 
the divergence in the kinetic energy of a single vortex in 
a two-dimensional Bose-Einstein condensate. However, 
the energy of a pair consisting of a meron together with 
another meron of opposite vorticity, or antimeron, is fi- 
nite since the topological deformation of the spin texture 
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FIG. 3: Merons in two dimensions, (a) A pair of merons with the same Pontryagin index but opposite vorticities. The arrows 
indicate the staggered magnetization, the periodic optical lattice potential is shown beneath and the core regions of the merons 
are shown as yellow discs, (b) The radius of the core region A is plotted as a function of temperature and imbalance. The core 
size diverges as the imbalance is reduced to zero since full rotational symmetry is then restored and merons no longer exist. 
The surface terminates at the critical temperature where the system becomes paramagnetic. 



cancels far away from the centre of the pair where the 
Neel vector always points in the same direction, as il- 
lustrated in Fig. E^a). Such meron-antimeron pairs are 
thermally excited below the critical temperature and are 
responsible for a Kosterlitz-Thoulcss phase transition as- 
sociated with the unbinding of the pairs [l^. Indeed, 
the entropic contribution of the merons above the criti- 
cal temperature is such that the system can lower its free 
energy through the proliferation of single merons. Topo- 
logical excitations are notably absent from the mean-field 
analysis and as a result our mean-field analysis becomes 
even qualitatively incorrect in two dimensions. To cor- 
rect for this, we used Monte Carlo results for a simi- 
lar system [19[ to estimate the critical temperature for 
the Kosterlitz-Thouless transition from the canted anti- 
ferromagnet to the paramagnetic phase. We found that 
this transition occurs at a significantly lower temperature 
than predicted by mean-field theory in two dimensions, 
as shown in Fig. 3) In particular, the phase transition 
occurs at zero temperature for the balanced case where 
the system rotationally symmetric. In three dimensions, 
however, we do not expect topological excitations or spin 
waves to significantly alter the phase diagram. Indeed, 
numerical studies yield only a 36% downward correction 
of the critical temperature in the balanced case [2^ . 

Outlook and conclusions — The imbalanced antiferro- 
magnet can be readily investigated experimentally with 
modern techniques. Imbalance is achieved in an ultracold 
Fermi gas by driving spin transitions with an RF field, 
and a Neel-ordered state can be realised by, for exam- 
ple, adiabatically ramping up the optical lattice poten- 
tial The Neel character can then be probed by the 
measurement of correlations in atom shot noise 21, 22| 
or using Bragg reflection 2^ 24 1 which also enables spin 
waves to be examined. The unbinding of superfluid vor- 
tices in two dimensions has been recently observed in 
a Bose gas [25| . and a similar interference experiment 
could in principle be performed to detect the Kosterlitz- 
Thouless phase transition here. Moreover, in situ imag- 
ing techniques based on scanning electron microscopy 
could be applied to explore topological and magnon ex- 
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FIG. 4: Phase diagram of the imbalanced antiferromagnet in 
two dimensions. The critical temperature for the canted phase 
(C) is drastically reduced compared to the mean-field result 
(dashed line) due to fluctuations and becomes a Kosterlitz- 
Thouless phase transition. Inset: a close-up of the Kosterlitz- 
Thouless transition for small imbalance. 



citations by resolving single lattice sites 2^ 27|. Finally, 
merons possess an internal Ising degree of freedom associ- 
ated with their Pontryagin index that is very suggestive 
of their application as a qubit in topological quantum 
computation, and this is a particularly interesting topic 
for further investigation. 

Methods — At half filling and for fc^T < U, with 
U being the energy cost of adding a second particle to 
a singly-occupied lattice site, the system is in the Mott- 
insulator phase and its low-lying excitations are described 
by the effective antiferromagnetic Heisenberg Hamilto- 
nian. Imbalance in the number of atoms of each spin 
can be incorporated by means of a constraint, which en- 
forces that the average total spin is equal to the desired 
magnetisation. Thus, the Hamiltonian of the imbalanced 
system is 77 = {J/'2)J2{i,j) Si-Sj-X^^ B-(Si-m), where 
hSi is the spin-i operator on site i, {i,j) denotes a sum 
over nearest-neighbour sites of the bipartite hypercubic 
lattice of interest here, and the effective magnetic field 
B acts as a Lagrange multiplier. The exchange constant 
J arises from the superexchange mechanism and is posi- 
tive. That is, the system can lower its energy by virtual 
nearest-neighbour hops only when there is antiferromag- 
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netic ordering. Note that the total spin operator ^ ■ 
commutes with H and thus the magnetisation is a con- 
stant of the motion. 

Within the usual mean-field analysis we obtain the free 
energy per site 
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where B^ (^b) = B — Jzm ± Jzn and we take the up- 
per (lower) sign for the A (B) sublattice, T is the tem- 
perature, kB is the Boltzmann constant, and n is the 
staggered magnetisation which is related to the aver- 
age on-site value of the spin by {Sa{b)) = m ± n. 
The Lagrange multiplier B is then found from the con- 
straint 9/(n, m; B)/9B ~ 0. Minimising the free en- 
ergy subject to this constraint, we find the critical tem- 
perature at which (n) becomes nonzero to be Tc — 
Jzmz/[2kB a.rcta.nh{2mz)]- At this point and through- 
out the antiferromagnetic phase, the constraint is solved 
by B = 2 Jzm and (n) is perpendicular to m in the min- 
imum of the free energy. By contrast, in the Ising model 
where n and m are restricted to be parallel, the critical 
temperature one obtains is Tc = zJ{l/A — m1)/kB- 

We perform a variational calculation substituting the 
spin-texture of the meron into Eq. ([T]) to determine 
its characteristic size for a range of temperatures and 
imbalance. A meron spin texture in the x-y plane 
with vorticity Uy — ±1 can be described by n = 
{ i/n^ — [nz(r)]2 cos 0, Uy \J t? — \nz(r)Y sin 0, (r ) } , 
where n = |(n)| and is the azimuth in the x-y plane. 
Exact meron solutions that follow from minimising the 
non-linear sigma model action with this texture have an 
out-of-plane component that behaves as n — n^(r) oc 
near the origin and that decays exponentially far from 
the origin [28|. A suitable variational ansatz that has 
the right properties near the origin but decays only 
as a power law for large r is riz(r) = n/[(r/A)^ -I- 1]^. 
In this case, the energy of a single meron due to 
the gradient term in the non-linear sigma model is 
Jn27r[(511/60) - tt - 61n2 + (1/2) ln(A/7rA2)] + 0[A-^) 
which diverges as the logarithm of the area A of the 
system. Incorporating the correct exponential decay 
for large distances changes the constant terms in this 
expression but leaves the logarithm unaffected, indicat- 
ing that our ansatz is indeed suitable to determine the 
characteristic core radius A. Minimizing the sum of the 
gradient energy and the integrated on-site free energy in 
Eq. ([3]) of the meron texture leads to the optimal size 
shown in Fig. [3ljb). The Pontryagin index here is +1/2 
and a texture with index —1/2 produced by substituting 
Tizir) — > —nzir) will produce an identical variational 
result, as noted earlier. Also, the core structure will play 
no role in the interaction so long as the spacing of a pair 



of merons is larger than A. 

The Kosterlitz-Thouless critical temperature has been 
calculated for an anisotropic 0(3) model using Monte 
Carlo techniques 19]. By making an analogy between the 
free energy contribution of our non-linear sigma model 
and the anisotropy in the anisotropic 0(3) model, we can 
use the numerical results of the anisotropic 0(3) model to 
estimate the critical temperature for the imbalanced anti- 
ferromagnet that incorporates fluctuations beyond mean- 
field theory. Although for rtiz ^ 0.2 our anisotropy is 
more complicated than that used in the Monte Carlo 
simulations, we do not expect a rigourous calculation to 
deviate significantly from the results in Fig. (4] 
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